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Abstract 

Dynamic nuclear spin polarization (DNP) mediated by paramagnetic point defects in semiconductors is 
a key resource for both initializing nuclear quantum memories and producing nuclear hyperpolarization. 
DNP is therefore an important process in the field of quantum-information processing, sensitivity-enhanced 
nuclear magnetic resonance, and nuclear-spin-based spintronics. DNP based on optical pumping of point 
defects has been demonstrated by using the electron spin of nitrogen-vacancy (NV) center in diamond, and 
more recently, by using divacancy and related defect spins in hexagonal silicon carbide (SiC). Here, we 
describe a general model for these optical DNP processes that allows the effects of many microscopic pro¬ 
cesses to be integrated. Applying this theory, we gain a deeper insight into dynamic nuclear spin polarization 
and the physics of diamond and SiC defects. Our results are in good agreement with experimental obser¬ 
vations and provide a detailed and unified understanding. In particular, our findings show that the defects 
electron spin coherence times and excited state lifetimes are crucial factors in the entire DNP process. 
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I. INTRODUCTION 


Point defects in solids are promising implementations of quantum bits for quantum comput- 
in^IEl. In particular, the negatively charged nitrogen-vacancy defect (NV center) in diamoncP has 
become a leading system in solid-state quantum-information processing because of its unique mag¬ 
netooptical properties, including long spin coherence timed^and the ease of optical initialization 
and readout of its spin state,l^even non-destructivel)^^. Since it has a high-spin electronic structure 
similar to the NV center in diamond, the divacancy in silicon carbide (SiC) has also been proposed 
to serve as a solid-state quantum biP. In fact, SiC hosts many other color centers that may also 
act as quantum bit^mni 

Recent demonstrations have shown coherent manipulation of divacancy 
and related defect spins in 4H-EI, I and SC-SiCP^. Coherent control of the electronic spin 

of the negatively charged Si vacancy has also been investigatecP^^i^^. Further milestones on the 
path towards robust SiC-based quantum-information technology have been the findings that iso¬ 
lated divacancy qubits have ~1 ms at low temperature^ and that isolated Si-vacancy qubits can 
operate at room temperatureP'. 

Coherent control of the electron spins of paramagnetic point defects makes it possible to 
control and manipulate other spins in the vicinity of the point defect. For instance, the prox¬ 
imate nuclear spins of the NV center in diamond can be polarizecP^, which can be a basis 
for quantum memoried^^®^, entanglement-based metrological device^, and solid-state nuclear 
gyroscoped^^I^. A recent demonstration has shown that nuclear spins proximate to divacancies 
and related defects in 4H- and 6H-SiC can be effectively polarizecP, an important step towards 
enabling long-lived quantum-information processing in this technologically mature semiconductor 
material. The transfer of the point defects’ electron spin polarization can also lead to hyperpolar¬ 
ization of the host material, thereby enabling sensitivity-enhanced nuclear magnetic resonance and 

spintronic applicationP^®o| 

Many of these applications rely on dynamic nuclear spin polarization (DNP) to mediate polar¬ 
ization transfer from the electron spin to neighboring nuclear spins through the hyperfine inter¬ 
action. Therefore a fundamental understanding of DNP processes is an important aspect of the 
technological development of nuclear spintronics. 

Jacques et developed an insightful spin Hamiltonian model that describes the polarization 
process of ^^N nuclei at the avoided crossing of the diamond NV centers spin sublevels, otherwise 
known as the level anti-crossing (LAC). The transverse part of the hyperfine interaction is respon- 
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sible for the exchange of electron and nuclear polarizatiorP^, while the spin selective non-radiative 
decay of the electron from the excited state (ES) is responsible for the maintenance of the electron- 
spin polarization. Continuous cycling of optical electronic excitation, flip-flops of electronic and 
nuclear spins, and non-radiative decay results in a polarization of both the electron and the nuclear 
spin populations. The so-called “excited-state level anticrossing (ESLAC) mechanism” is when 
spin flip-flops predominantly occur in the excited state. At larger magnetic fields, when the EAC 
occurs in the ground state (GS), the analogous ground-state level anticrossing (GSEAC) process 
occur^^ll22l Xo understand the particular features observed in the GS dynamic nuclear spin polar¬ 
ization of a nuclei adjacent to the vacancy of the NV-center, Wang et a/.^ recently developed 
a model capable of describing general hyperfine interactions, i.e. for nuclei not on the symme¬ 
try axis. Eurthermore, to take into account the effect of external strain and spin relaxation and 
dephasing processes, Eischer et recently proposed a density-matrix model. 

While these models capture certain phenomena of DNP, a more general model showing pre¬ 
dictive power over several color-center systems would be an important development. Here, we 
propose an extended model that (i) handles general hyperfine interactions of paramagnetic defects 
with symmetrically or non-symmetrically placed nuclei spins, (ii) takes into account simultaneous 
GSEAC and ESEAC processes, and (iii) tracks the evolution of spins with time explicitly param¬ 
eterized. Our model provide insight into the phenomena of the dynamic nuclear spin polarization 
mechanism for both the NV center in diamond and the divacancies in SiC, and explains several 
experimental observations. Throughout our investigation, electron-spin decoherence appears as 
an important limiting factor for the polarizability of the nuclear spins. We show that considering 
electron-spin coherence will be vital to maximizing the performance of DNP in practical applica¬ 
tions. 

In Section |I^ we briefly describe the electronic structure and the corresponding spin properties 
of the ground and excited states of the considered point defects, namely, the NV center in diamond 
and the divacancy defects in 4H- and 6H-SiC. These defects’ structures are then used to define the 
spin Hamiltonian, which is given together with the model of the dynamic nuclear spin polarization 


in Section In Section [1^ we also summarize the basic concept and parameters of the model 

we 


to calculate the nuclear spin polarization as a function of different variables. In Section IV 


describe the ab initio methods and models to calculate the spin related properties of NV center in 
diamond and divacancy defects in SiC, along with the corresponding results. The full hyperfine 
tensors of the studied nuclear spins are calculated both in the ground and excited states and are 
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important parameters in the DNP model. In Seetion |V| we provide the results and an analysis of 
DNP for the NV eenter in diamond and the divaeaney in SiC. Finally, we summarize our findings 
in Seetion rvTl 


II. ELECTRONIC STRUCTURE OE NV CENTER IN DIAMOND AND DIVACANCY IN SIC 

The geometry and the eleetronie strueture of the negatively eharged NV eenter in diamond 
have been diseussed previously based on highly eonvergent ab initio plane wave large supereell 
ealeulationd^SmSl diamond NV eenter is a eomplex that eonsists of a substitutional nitrogen 
adjaeent to a vaeaney in diamond and possesses symmetry (see Fig. [^. The defeet exhibits 
a fully oeeupied lower oi level, and a double-degenerate upper e level filled by two parallel-spin 
eleetrons in the gap with an S' = 1 high-spin ground state. The S' = 1 exeited state is well 
deseribed by the promotion of an eleetron from the lower defeet level to the upper level in the 
gappl. The eleetron spin may interaet with nuelear spins: ^^N or ^®N possessing J=lorJ = l/2 
nuelear spin, respeetively, or with I = 1/2. 

The hyperfine interaetion between the eleetron spin and nuelear spin has been studied by means 
of ab initio methods in previous publieations, both in the and in the E;^. In the 

GS, the eleetronie spin-spin dipole interaetion eauses the fine eleetron strueture to have a zero- 
field splitting, where “zero field” refers to zero external magnetie field. In high purity and low 
strain diamond samples, this splitting ean be deseribed by a single parameter, Dqs=2.S1 GHz, 
whieh separates the m 5 = 0 and the ms = ±1 sublevels within the S' = 1 manifold. At room 
temperature, the fine strueture in the eleetronie ES shows a similar feature, exeept that its zero- 
field splitting is only Z1 es=1-42 GHz. At lower temperatures the fine strueture beeomes more 
eomplieated, hindering off-resonantly pumped DNP in the EI^. Between the ES and GS triplets, 
non-radiative relaxation pathways through singlet states seleetively flip ms = ±1 states to ms = 
0 state in the optieal exeitation eyele!^^®, whieh allows optieal spin-polarization of NV eenter 
(Eig.Q. 

We also study the divaeaney defeets in SiC. Before deseribing the eleetronie strueture of the 
divaeaney, we will briefly diseuss the host semieonduetor. SiC has about 250 known different 
polytypes, whieh share the same basal hexagonal lattiee but have different staeking sequenees of 
Si-C bilayers perpendieular to this plane. The most industrially important polytypes are the 4H 
and 6H poly types. The inequivaleney of the erystal planes leads to the so-ealled h and k types 
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FIG. 1. (Color on-line) A schematic diagram of (a) the structure showing the Cs rotation axis and (b) the 
electron configuration of the NV center in diamond, (c) The magnetic field dependence of fhe ES and GS 
spin-sublevel energies, showing fhe ESLAC and GSEAC. (d) The NV cenfer’s opfical polarizafion cycle 
and (e) fhe nuclear spin polarizafion cycle. In (c) and (d), L>gs=2-87 GHz, Oes=E42 GHz are fhe zero-field 
consfanfs. In (d), fhe green dashed arrows represenf fhe non-radiafive decays, which are mediafed by spin- 
orbif couplings and vibrations. The fhick grey arrows represenf fhe opfical absorpfion/emission pafhs. The 
wavy lines represenf phofon absorpfion/emission in fhe visible (red) and near-infrared (brown) regions, (e) 
Af zero magnefic field {B = 0), fhe |0 i) level is separafed by fhe zero-field consfanf from fhe | — 1 t) level. 
Applying a B = Hlac > 0 field causes fhe fwo sfafes fo form an avoided crossing, where fhe small gap 
is infroduced by fhe hyperfine inferacfion (A). In fhis condifion, hyperfine coupling (brown circular arrow) 
and non-radiafive decay (green doffed line) are responsible for fhe nuclear spin polarizafion in fhe opfical 
cycle. 
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of bilayers in 4H SiC, and h, ki, ^2 types of bilayers in 6H SiC. In the inequivalent bilayers, the 
erystalline environment in the seeond, third, ete. neighborhood will be different, yielding a sim¬ 
ilar but quantitatively distinguishable eleetron strueture for the eorresponding point defeets. For 
divaeaney defeets, adjaeent silicon and carbon atoms are absent. In 4H SiC, the four inequivalent 
forms of divacancy are the hh, kk, hk and kh configurations. The first two (axial) configurations 
have symmetry, and the second two (basal-plane oriented) configurations have Cih symmetry. 
We focus our study here on the axial configurations, whose symmetry is the same as that of 
the NV center in diamond. In 6H-SiC, there are three axial configurations of divacancies: the 
hh, kiki and ^ 2^2 configurations. In 4H-SiC, the hh and kk configurations have been already 
associated with zero-phonon-line photoluminescence peaks and spin transitions, measured with 
optically detected magnetic resonance (ODMR) signals and sometimes called the PLl and PL2 
center^2ISI113 ^n additional center has been found, PL6, which also exhibits ODMR and similar 
physical properties to the axial divacanciej '^ l ^^*^^l The physical structure of PL6 has not yet been 
identified. In 6H-SiC, the PL and ODMR lines that have been labeled QLl, QL2, and QLd^are 
also associated with the axial divacancies (Section [TV]), namely the kiki, hh and ^ 2^2 configu¬ 
rations, respectively. All these axial divacancies share the same electronic structure depicted in 
Fig.|^ The carbon dangling bonds of the Si-vacancy part of the defect create a double-degenerate 
e-level close the valence band edge, which is occupied by two electrons with parallel spins. Thus, 
the neutral divacancy has a high spin (S' = 1) ground state. The excited state may be described as 
the promotion of an electron from the lower defect oi level to this e leveP^, akin to that of the 
NV center in diamond. In SiC crystals, beside isotopes (with a 1.3% natural abundance), ^®Si 
isotopes (with a 4.7% natural abundance) have 1 = 1/2 nuclear spins that may interact with the 
divacancies’ S' = 1 electron spins. Very little is known about the nature of the triplet ES and the 
dark singlet states responsible for the electron spin polarization for divacancy in SiC. However, 
recent measurements indicate!^ that the triplet ES has a similar electronic structure to that of NV 
center in diamond. Experimental results imply symmetry in the ES of axial divacancies even 
at low temperature^. In this work, we assume that the SiC divacancy has a similar model for the 
optical spin polarization processes as does the NV center in diamond. The measured Hgs and Des 
zero-field splittings of axial divacancies in 4H and 6H polytypes are summarized in Table 
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FIG. 2. A schematic diagram of the structure of the divacancy in SiC. The defect levels in the gap as well 
as the corresponding ground and excited states are shown. The high energy upper empty e level does not 
play a role in the excitation process. Those Si atoms that participate in DNP are labeled by Sina and Sinb. 
The similarity between the electronic structure of the NV center in diamond (c.f., Fig.[T]) and that of the SiC 
divacancy is apparent. 

III. MODEL OF THE DYNAMIC NUCLEAR SPIN POLARIZATION PROCESS 

A. Modeling the dynamic nuclear spin polarization in the optical cycle 

The degree of nuelear spin polarization (P) is the main observable in DNP measurements. Un¬ 
derstanding its dependenee on external variables like magnetie field and temperature, and internal 
variables like the details of the hyperfine tensor, has great importanee from an applieations point 
of view. A theoretieal model is a erueial eomponent of this understanding. 

In the fortheoming seetion, we briefly review the model of Jaeques et al^on the dynamie nu¬ 
elear spin polarization of the NV center’s nitrogen nucleus and extend it by a generalized derivation 
of the basic equations. Next, we describe our model that takes into account many microscopic fea¬ 
tures and processes that have been overlooked or not integrated in previous models. The most 
important features that we integrate are the electron spin decoherence in the ground and excited 
state, the short lifetime of the excited state, the anisotropy of the hyperfine tensors, the angle of the 
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external magnetie field, and the overlap of the ground and exeited state’s spin flipping proeesses. 


1. Previous model on dynamic nuclear spin polarization 


The dynamie nuelear spin polarization proeesses are sensitive to the details of the spin Hamilto¬ 
nian of the eonsidered point defeet and the optieal eleetron spin polarization eyele. In the simplest 
ease, e.g. an NV eenter in diamond and an adjaeent 1 = 1/2 nuelear spin loeate on the symmetry 
axis of the defeet, the proeess ean be understood as a two-step meehanism. Continuous optieal 
exeitation polarizes the eleetron spin in the Ms = 0 spin state, while the nuelear spin may be any 
linear eombination of the spin-up |t) and spin-down IJ,) states. Near the vieinity of a LAC, the 
hyperfine interaetion eouples the nuelear and eleetron spins effeetively and rotates the two spin 
state |0 J,) into the state | —11)- This state is then transformed into the state |0 t) by the optieal 
exeitation and deeay proeesses shown in Fig.[^(e). As a result, the nuelear spin eomponent IJ,) is 
flipped and the |t) is predominantly populated. 

The spin Hamiltonian that governs the nuelear spin flips is well known for the NV-eenter in 
diamond^. For a single nuelear spin adjaeent to the defeet, the ground-state Hamiltonian of the 
eleetron-nuelear spin system ean be written as 

Hgs = S^DgsS + /reB^geS -f S^Aqs! + (1) 


where S and I are the eleetron and nuelear spin operators, Dqs and Aqs are the tensors of zero-field 
and hyperfine interaetion in the ground state eleetron spin eonfiguration of the defeet, respeetively, 
B is the external magnetie field, ge and gN are the g tensors of the eleetron and nuelear spin, and 
Pb and /xn are the Bohr and nuelear magnetons, respeetively. For the eonsidered defeet, the g 
tensors are nearly isotropie and simplify to a sealar value, sueh as g^. = 2.0023. 

For the NV eenter, the zero-field interaetion term the first term on left hand side of Eq. ([^, 
ean be written as 

4fs = S^DgsS = Dos - 

The hyperfine interaetion term Hhyp, whieh is the third term on the left hand side of Eq. ([^, 
eouples the eleetron and nuelear spins and therefore plays a key role in DNP. When the nuelear 
spin resides on the symmetry axis of the defeet the hyperfine tensor A is diagonal with diagonal 
elements A± and Ay. The hyperfine interaetion term is written ad^ 




AgsI 


A 


GS 


SM. + S.D 




(3) 
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where S± and I± are the eleetron and nuelear spin ladder operators, respeetively, and Sz and are 
the 2 ; eomponents of the eleetron and nuelear spins, respeetively. The first term on the left hand 
side of Eq. Q is responsible for the flipping of the nuelear and eleetron spins governed by A±. 

The NV-eenter in diamond possesses symmetry in its exeited state (ES) at elevated 
temperature^^, and therefore the exeited state Hamiltonian has the same form as Eq. ([^. How¬ 
ever, the zero-field-splitting tensor Des and the hyperfine tensor Aes differ from those in the 
ground-state eleetronie eonfiguration. 

Eor the sake of a general deseription, we utilize the density matrix formalismP^ to derive the 
steady state nuelear spin polarization of the dynamieal proeess. Eirst, we express the wave funetion 
of an eleetron and nuelear spin system in a basis. 


^ = '^Cmn(,t)\mn), ( 4 ) 

m,n 

where \mn) = |m) ® |n), m and n are the eleetron and nuelear spin projeetions on the quantization 
axis and Cmn are eoeffieients. The spin-density matrix of the system ean be obtained from the 
eoeffleients, 

while the nuelear spin-density matrix ean be obtained from the partial traee of Pmn,m'n', 

*hnn'(f) ~ ^ Pmn,mn'{t') ■ (6) 

m 

To deseribe the time evolution of the diagonal elements of the nuelear spin density matrix in the 
dynamieal proeess, for timeseales larger than the average length of a dynamieal eyele, we write 
the kinetie equations in the form 


^nn = C{J) 


(7) 


where the dot represents time differentiation. The seeond term on the right hand side deseribes 
the nuelear spin relaxation due the environment. Note that all the nuelear spin projeetions are 
assumed to relax equally with rate p. The first term on the right hand side deseribes the ehange of 
the nuelear spin projeetion due to the interaetion with the eleetron spin and the external magnetie 
field. c(J) is the number of optieal eyele per unit time and is the averaged variation of 

the nuelear spin projeetion during the free evolution time, i.e. between two optieal exeitations. 
This last term ean be determined from the spin Hamiltonian of the system. Note that both terms 
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depend on the intensity J of the exeitation laser. At low intensities, the rate of optical excitation 
c( J) depends linearly on the intensity. For strong laser excitation, when the rate of excitation may 
be comparable with the periodicity of the spin rotations due to hyperfine interaction or transverse 
magnetic field, the time average of depends on J too. In the following, we assume weak 
laser intensities, and thus only c (J) depends on the intensity. Generally, A^nn can be written as 

A^nn ^ Pn'n^n'n' ^nn^^Pnn '7 ( 8 ) 

n' n' 

where pnn' is the average probability of flipping the nuclear spin from \n) to \n'). The first and 
second terms on the right hand side describe the probabilities of flipping the spin in and out of 
the state |n), respectively. We note that for the Hamiltonian specified in Eq. 0-0, n' = n±l 
hold^. In the following, we restrict our derivation to the case of / = 1/2, however, we do not 
impose any other losses in the generality of the spin Hamiltonian. In this case, Eq. 0 reads as 






■Pt^± 


-±- ’ 
2-^2 


(9) 


where p± are the average probabilities of raising and lowering of the nuclear spin projection be¬ 
tween two optical cycles. 

When dynamic nuclear spin polarization is in a stationary state, the different spin rotation 
processes are balanced and: 

<i), 1 , 1 = <i)_i_i. (10) 

2 2 2 2 

By using the definition P = $ , i , i — <h_i_i and the normalization condition , i , i + i_i = 

■'■ 2"'"2 2 2 ^ 2"''2 2 2 

1, from Eq. ([T0|) we can express the steady state nuclear spin polarization as 


P+-P- 

P+ + P- + K 


( 11 ) 


where n = rijc (J). 

As can be seen, the flipping probabilities and p_ play an important role in this process. In 
the previous modeP^, these were determined analytically from the simplified spin Hamiltonian 
Eq. 0-0 for the case S = 1 and / = 1/2, as 


p^{B) = 2\{0i\+)\l\{-lt\+)\l, (12) 

P-{B) = p+{-B) , (13) 

where |-|-) is the eigenstate of HEsiB). This state is a mixture of |0 j,) and |—11) states due to 
hyperfine coupling. 
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2 . Our extended model 


We now discuss the details of our extended model, which is generally applicable and shows 
predictive power. First, we apply modifications to the spin Hamiltonian in Eq. Q. 

In previous models, the external magnetic field was aligned parallel to the axis of the de¬ 
fect. However, experiments show that DNP strongly depends on the misalignment of the magnetic 
fielcP^. To describe the effect of this misalignment, we allow small deviations of the direction of 
the magnetic field from the symmetry axis of the defects. Furthermore, the nuclear Zeeman effect 
was also neglected in previous models. At stronger magnetic fields, such as at 500-1000 Gauss, 
the nuclear Zeeman splitting can reach few MHz, which can be comparable to the hyperfine inter¬ 
action. We thus take this effect into account in our model. 

In the general case, the nucleus with non-zero spin is not on the symmetry axis of the defect. In 
this case, the symmetry of the spin Hamiltonian is reduced, i.e. the hyperfine tensor A may have 
three non-degenerate eigenvalues A^x, Ayy, and A^z, and the eigenvector a^, which corresponds 
to the eigenvalue Azz, may have a non-zero angle of 9 with the symmetry axis. The azimuthal 
angle cp may be chosen to zero without limiting the generality. The effects of symmetry-breaking 
hyperfine interactions have been included in previous considerations to some extenlP^*^. However, 
a consistent description of the DNP process with a general hyperfine tensor has not been carried 
out so far. 

In our model, we consider non-diagonal hyperfine tensors, which can be parameterized by their 
eigenvalues and angle 9 as follows: 

iThyp = S^Ai = (us) Adiag(ui) , (14) 


where U describes a rotation that transforms the Cartesian basis to the eigenbasis of tensor A, 
and Adiag = UAU^ is the diagonal tensor of elements Axx, Ayy, and Azz- Note that for a general 
hyperfine tensor the spin Hamiltonian may contain S±Iz, Szl±, and S±I± terms that allow a wide 
range of spin rotation processes to occur (see Appendix for spin Hamiltonian matrices). 

The nuclear spin flipping probabilities play a key role in the determination of the polarization. 


defined in Eq. 11 The main innovation of our model is that it uses different definitions for the 
probabilities p+ and p_ than previous model^ and includes important effects from the excited 
and ground states’ spin Hamiltonians as well as external driving forces from the surrounding spin 
bath. In the rest of this section, we describe the main concept of our new considerations. 
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For simplicity, in this section we restrict ourselves to the case of positive external magnetic 
field, whieh shifts the energies of Ms = +1 and M 5 = — 1 levels upward and downward, respee- 
tively. In sueh a case, at vieinity of -Blac> the Ms = —1 spin state mixes with the Ms = 0 state. 
For simplieity, we do not eonsider the interaetion of the Ms = +1 state with other states here. 
However, we include it in our later calculations (see Seetion [riIB| ). 

Nuelear spin rotation ean occur both in the electronic ground and excited states. The interplay 
of these rotations and non-radiative, spin-seleetive electronic decay is responsible for the nuclear 
spin polarization. In the most general ease, to aehieve a net driving force toward nuelear spin 
polarization, the ground-state and the subsequent excited-state spin rotation proeesses have to 
fulfill the following eriterion: starting from the Ms = 0 eleetron spin state, the electron spin may 
be flipped into the Ms = —1 state, but the starting nuclear spin state must be flipped into the 
opposite spin projection. When electron spin flip-flops occur, the M5 = — 1 electron spin is then 
transported into the initial M 5 = 0 state by non-radiative decay. Optical cycles ean therefore 
flip nuclear spins. When the rates of the flipping processes |0 j,) — )■ |0 t) and |0 t) — )■ |0 j,) are 
different due to some sort of asymmetry, the repetition of these dynamical cycles induces different 
population of the nuclear spin states and non-zero polarization. 

All the possible spin rotation processes that fulfill the aforementioned eriterion and are allowed 
by the general spin Hamiltonian (see Appendix for details) are depicted on a sehematie model of 
the dynamical cycle in Fig. The probability of the ground or excited state spin rotation pro- 
eesses p^‘(xinitiai| XKnai), where Xinitiai and XKnai represents the initial and final spin configurations, 
respectively, ean be determined by the spin Hamiltonian of the defeet in the eonsidered states ’St’. 
The probability of spin flip in a joint ground state - excited state spin rotation proeess is then the 
product of the probabilities of the two separate rotations. With the above requirement, the prod¬ 
ucts p°^( 0 , ± 1 / 2 | Xlnter)P^^(Xlnter| “ 1 , =F 1 / 2 ) and p^^( 0 , ± 1 / 2 | Xinter) ( Xinter | 0 , =F 1 / 2 ) define 

the probability of spin rotations that induee driving forces toward nuclear spin polarization, where 
Xinter represents an intermediate spin eonfiguration and |±1/2) is a eomprehensive notation for It) 
or It) nuelear spin states. 

In eontrast to previous models, our two-state evaluation model requires suitable spin rotation 
mechanism for DNP not separately but simultaneously in the ground and excited states. For ex¬ 
ample, at i?GSLAC 7 the ground-state hyperfine eoupling flips the nuelear and eleetron spins with 
high probability. This proeess is represented by orange arrow on the top panels of Fig. Optical 
excitation transports this spin state to the exeited state, where it starts to evolve in aeeordanee with 
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FIG. 3. (Color online) A schematic diagram of the evolution of a coupled electron and nuclear spin system 
through the four steps of a complete dynamic nuclear spin polarization cycle. The depicted spin rotation 
processes, shown by different paths along the black arrows, preserve the spin projection Ms = 0 of the 
electron spin but flip the nuclear spin (a) from down to up, and (b) from up to down, in a four-step mech¬ 
anism. The states of up and down nuclear spin projections are depicted with blue and green backgrounds, 
respectively. In the general case, the interplay of all of these parallel processes determines the net polariza¬ 
tion of the nuclear spins, with the process rates defined by the ground and excited state spin Hamiltonians 
(see text and Appendix for more information). The upper and lower path of arrows with orange (thick 
light grey) outline show the most prominent spin rotation processes at the GSLAC and ESLAC for positive 
magnetic fields, respectively. The background spin bath, which works against DNP, is represented by grey 
background. 


the new spin Hamiltonian of the exeited state. To obtain a high probability for suitable eleetron 
and nuelear spin flips from the two sueeessive time evaluations, the spin state must be unehanged 
in the exeited state. For general hyperfine tensors of the exeited and ground state that are not nee- 
essarily identieal, this eondition may not be the ease. The role of minority proeesses, represented 
with blaek arrows in Fig. inereases as the symmetry of the system beeomes more distorted. 
Examples of this sort of distortion inelude a misaligned magnetie field or hyperfine interaetions 
with low symmetry. Sueh effeets ean lower the polarizability of the system or eause unexpeeted 
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resonance effects, which can be captured by our two-state evaluation model. 

The next new consideration in our model is the determination of the preferential direction of 
the nuclear spin polarization. The DNP process is in a stationary state when the polarization and 
depolarization processes are in equilibrium (i.e. Eq. ( [T0| ) is satisfied). This state corresponds to 
a certain electron and nuclear spin configuration, which is not necessarily equal to one of the 
energy eigenstates of the ground or excited-state spin Hamiltonian. Generally, the nuclear spin is 
polarized or preserved in a state that has the longest lifetime in the dynamical cycle. When the 
hyperfine interaction is isotropic or symmetric, the energy eigenstate |0 t) has the longest, 
quasi-infinite lifetime, since in this state there is no hyperfine coupling between the electronic 
and nuclear subsystem^. On the other hand, for general hyperfine interactions of the excited and 
ground state, the longest-lived spin state, which is the stationary spin state, can be a non-symmetric 
stateP. 

This non-symmetric state can be written as a linear combination |0) 0 (a |t) + /9 |i))- If the 
nuclear spin state is the eigenstate of the nuclear spin operator Je of quantization axis e, then unit 
vector e of the three dimensional space shows the preferential direction for nuclear spin polariza¬ 
tion in the dynamical cycle. In our model, we determine this e direction for every different set of 
parameters. To find this direction, we used the following condition: the nuclear spin state Ite)^ 
which is the eigenstate of Je with eigenvalue -1-1/2, has the longest lifetime when its precession is 
minimal. The precession would result in a continuous transition between the two eigenstates of Jg, 
i.e. Ite) lie), which then reduces the lifetime of the nuclear spin state and lowers polarization. 
Since the system spends most of the time in the ground state during the dynamical cycle, we as¬ 
sumed that the ground-state precession should be minimal in the stationary state of the dynamical 
cycle. This means that we minimized 0, teI 0, ie) with respect to direction e. At the minimum, 
we find the direction where the nuclear spin state |te) is preserved for the longest time. This state 
is therefore the stationary nuclear spin state in the dynamical cycle. 

In the optical cycle, the length of the free evolution of the spin system is limited by the lifetime 
and decoherence of the electron states (see below). When the lifetime of the ground or excited state 
is shorter than the characteristic time of the spin rotation processes, the oscillatory probability is 
not averaged out. In extreme cases, the spin-flip probability is effectively reduced by the short 
lifetime of the electron state (see Fig. for the case of the NV center in diamond in Section fVA] ). 
This case is more pronounced in the excited state, since the excited state lifetime of both the 
NV center in diamond and the divacancy in SiC is only 10-15 ns. For a hyperfine interaction 
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of 10 MHz, the oscillation time of the nuclear spin-flip probability is on the order of 100 ns 
at -Blacj is much longer than the lifetime of the excited state. In contrast to previous models, 
we calculate the probabilities by explicitly taking time into account in the spins’ excited-state 


evolution (see Section |IIIB| for more detail). We assume that the evolution time in the excited 
state is exponentially decaying. The characteristic time of decay is connected to the lifetime of the 
excited state (tes) determined by experiment. On the other hand, we assume that the oscillatory 
probabilities are averaged out in the ground state, where the system spends sufficient time for this 
to happen. 

Finally, we include additional important effects in our model: electron spin dephasing and 
spin-lattice relaxation effects. The surrounding spin bath of nuclear spins and paramagnetic de¬ 
fects disturbs not only the nuclear but the electron spin of the considered spin systerrP^, causing 
decoherence of the superposition states. This decoherence can be described by the Lindblad equa¬ 
tion, 

dp 


dt 


i 

h 


H,p + Lp, 


(15) 


where p and H are the density operator and the Hamiltonian of the considered system, respectively. 


The final term on the right hand side of Eq. (151, which describes the effects of the environment, 
can be written as 

ip = {(?„(?*,p}), ( 16 ) 


where 7 „ are the decay rates of processes described by the Lindblad operators (7„. Decoupling 
of the superposition state of the |i) and \j) states can be taken into account by the operator C = 
|i)(z| — \j){j\. In our case, the state |- 1 -) = a |0 ),) - 1 - | —11), which is responsible for the most 

prominent nuclear spin flipping process, relaxes with T 2 characteristic time due to the electron 
spin dephasing processes can be accounted by the Lindblad operator |0)(0| — | —1)(—1|. Besides 
dephasing, electron spin-lattice relaxation can also hinder the nuclear spin flipping processes, by 
depopulating the hyperfine coupled states |0 1) and |—1 t). This effect can be described by the 
Lindblad operators \i){j\. Recent experiments have showiJ^that the electron spin coherence time 
T 2 in the excited state of the NV center in diamond is on the order of the excited-state lifetime. 
To take into account the aforementioned effects, we assume that the intact evolution time of the 
electron and nuclear spin system, described by the spin Hamiltonian in Eq. 0. is effectively 
reduced in the excited state due to the short coherence time of the electron spin. Therefore, we 
used a scaled excited state lifetime r^g = in our model, where z/ is a free parameter of 
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the property 0 < z/ < 1. By eonsidering eleetron spin deeoherenee and spin-lattiee relaxation 
elfeets, ~ 1 /^ 2 * + where Ti is the eharaeteristie time of the spin-lattiee relaxation. 

As dephasing is the fastest spin relaxation proeess, is assumed to be elose to, but somewhat 
smaller than T^. For the ground state, where the evaluation time is not ineluded explieitly, we 
sealed down the probability of nuelear spin rotation proeesses by a faetor of /i. 

Using the results of experiment and first-prineiples ealeulations for the parameters of the above 
deseribed model, three free parameters remain. All of them are related to some extent to the effeet 
of deeoherenee and spin-lattiee relaxation. Parameter k = r]/c{J) in Eq. ( [TT] ) is mainly related to 
the nuelear spin-lattiee relaxation, and the two parameters u and /i, introdueed above, are elosely 
eonneeted to the eleetron dephasing effeets in the exeited and ground state, respeetively. 

Despite these new eonsiderations, the model that we have deseribed still has limitations. Sinee 
it is a single eyele model, i.e. evolution of the spin system is taken into aeeount in a single optieal 
eyele, eomplieated proeesses that take plaee over many eyeles are not modeled. Sueh a meehanism 
appears for / > lEHi. 


B. Method of calculation of the nuclear spin polarization 


In this seetion, we speeify the equations that are used to ealeulate the nuelear polarization, P, 
in the framework of our model deseribed in the previous seetion. 


As Eq. (11) shows, the polarization ean be ealeulated from the probability of nuelear spin up 
and down flips (p+ and p_) in a dynamioal eyele, respeetively, and from the rate of spontaneous 
nuelear spin flips k, due to the baekground spin bath. The latter quantity is one of our model’s 
free parameters. Sinee the EAC happens for two values of the external magnetie field, ±i?LAC> the 
nuelear spin flipping probabilities ean be divided up into two parts. 


(~ 1) I /1 n\ 

P+=P+ +P+ , (17) 

(-1) , (-HI) 

p_ = p'_ ' ^ \ 


where and and and p^^^ represent the spin up and down flipping probabilities 
due spin rotation proeesses in the subspaee Ms = {0,-1} and Ms = {0, -1-1} , respeetively. Eor 
positive values of B the seeond terms on the right hand side of the equations have only minor 
eontribution. However, their role inereases as S —)■ 0. At 5 = 0, p+ and p_ beeome equal, and 
therefore, P\b=o = 0- our model, we use the approximatioEplp^^^(i?) = p’C^\—B) and 


16 





= p+ -B). The nuclear spin-flipping probabilities, corresponding to the Ms = {0, —1} 
subspace, can be defined as. 


p+ 11 Xi) -1 t)r + p^^(xi|o t)], (18) 

i 


p- tlxi) [p^^(xil -1 i)r + p^^(x*|o;)], 


where T is the probability of non-radiative decay from the electron spin state |±1) of the excited 
state to the ground state spin state |0) and Xi the final and initial states of the ground and 
excited states’ time evolution, respectively. The summation goes over states |0 t),|0 |), |—11)> 
and |—1 I), see Fig.|^ 

To evaluate Eq. (18), we define the probabilities p^'(xinitiai| XFinai)- For the ground state spin 
rotation processes, we averaged out | (xFinail x(f)) 1^ In time, as 


1 


f^GS 


P (Xinitiall XFinal) ‘^I.F rp / 

J-GS Jo 


(XFinal|e IXinitial) 


dt. 


(19) 


where we consider the integration time Tqs —)■ oo. This limit means that the system spends 
enough time in the ground state for the oscillatory probabilities to be completely averaged out. 
The factor afl takes into account the destructive effect of the electron spin decoherence and spin- 
lattice relaxation. is a two-value function that takes 1 when the initial and final states are the 
same and takes 0 < /i < 1 when spin rotation occur. The parameter p is a fitting parameter of our 
model. 

In contrast to that of the ground state, the excited state’s lifetime is short. In this case, the 
flipping probabilities strongly depend on the duration of the excited state’s evolution time, since 
it can be shorter than the periodicity of the oscillatory probabilities p^^{xi\ X/)> see Fig. For a 
proper description, we have to include time in our considerations. 

The flipping probabilities, correspond to the excited state’s time evolution, thus 


P (Xinitiall XFinal) / (XFinal|c 


Tes 




IXinitial) 


dt. 


( 20 ) 


Jo 

where g{t) is the probability distribution function of the effective length of the excited state’s 
evolution time. g{t) is assumed to be an exponential distribution 


g{t) 



■iA 


( 21 ) 
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where is the eharaeteristie time of the deeay. In our model, is the average effeetive time of 
evolution in the exeited state, whieh is eonsidered to be proportional to the exeited state’s lifetime 

"^ESj 

'^Es = ^'Tes- (22) 

Here, u, whieh is the last free parameter of our model, is a sealing faetor that takes into aeeount 
eleetron spin relaxation effeets that ean ehange the net evolution time. 

Finally, as diseussed during the deseription of the model, the preferential direetion e of the 
nuelear spin polarization may deviate from the axis of the defeeP^. We determine this direetion 
from the ground-state-spin Hamiltonian by finding the direetion e, where the nuelear spin state 
|te)> whieh is an eigenstate of Je with -1-1/2 eigenvalue, possess the lowest probability to flip into 
|J,e) in the ground state. This means that the preeession of state |te) is minimal and this state is 
therefore the stationary state of the dynamie nuelear spin polarization eyele. The eriterion used is 


minp^^(0 tel 0 te) • 


(23) 


The ealeulation of the nuelear spin polarization P for a given system defined by the spin Hamil¬ 
tonian of Eq. Q ean now be earried out by finding the suitable nuelear spin state |t) and |t) by 
Eq. ([23|), and using it to determine the flipping probabilities of the spin state in aeeordanee with 


Eqs. (19) and (20). Erom these values, the total probability of nuelear spin down-to-up and up-to- 


down flipping ean be obtained by Eq. (18) whieh provides the polarization of the nuelear spin via 


Eq. ([nj). 

The predefined parameters of the model are the diagonal elements and angles of the exeited 
and ground state’s hyperfine tensors, A^l, A^y, A^^, Oqs, A^l, A^y, Aff, and 6es, the zero-field¬ 
splitting parameter of the ground and exeited state’s zero-field-splitting tensors, Dqs and D^s, the 
rate of non-radiative deeay T, and the exeited state’s lifetime tes. The free parameters used to fit 
the theoretieal eurves to the experimental ones are /i, u, and k. 


IV. AB INITIO CALCULATIONS OF THE MODEL PARAMETERS FOR THE NV CENTER IN 
DIAMOND AND THE DIVACANCY IN 4H AND 6H SIC 

Sinee some key parameters of the DNP model have not been measured, we apply ab initio meth¬ 
ods to ealeulate them. In partieular, we ealeulate the full hyperfine tensor of seleeted proximate 
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J = 1/2 isotopes in the ground and exeited states. In addition, we identify the axial divaeaney eon- 
figurations in 6H-SiC by ealeulating their Dqs parameters in order to provide a direet eomparison 
for DNP proeesses in 6H-SiC. 

We earry out first-prineiples density-funetional theory (DFT) ealeulations to study NV eenter 
in diamond and axial divaeaneies in 4H and 6H SiC. We apply 512-atom supereell for diamond, 
and 576-atom and 432-atom supereells for 4H and 6H SiC, respeetively. We apply F-point sam¬ 
pling of the Brillouin-zone, whieh suffiees to ensure eonvergent eharge and spin densities. We 
utilize the plane wave basis set together with projeetor augmented wave method as implemented 
in VASP5.3.5 eodeP^. We apply our in-house eode to ealeulate the GS zero-field splitting from 
DFT wave funetion^^fi^, whieh has been tested and shown to provide good results with using 
Perdew-Burke-Emzerhof (PBE) funetionai^. We apply the HSE06 hybrid funetionaP^'^ to eal- 
eulate the hyperfine tensors of seleeted J = 1/2 nuelei where the spin polarization of the eore 
eleetrons are taken into aeeounP^. We apply the eonstrained DET method to ealeulate the ES spin 
densit)^^. 


A. Results on NV center in diamond 

The full hyperfine tensors of and eoupled to NV eenters in diamond have not been 
experimentally determined. We thus apply ab initio ealeulations to obtain these parameters. The 
GS-hyperfine tensors of NV eenter in diamond have been reeently eharaeterized in detaiP^. In 
partieular, we analyze the hyperfine tensors in the GS and ES for eoupled ^®N nuelei, for whieh 
the DNP was thoroughly studied experimentall)^^. We provided the ES hyperfine tensor by PBE 
funetionaP^, where we showed that the hyperfine eonstants are anisotropie for ^^N. We list the 
eorresponding HSE06 values in GS and ES in Table 


B. Results on divacancies in 4H and 6H SiC 

We ealeulate the eleetronie strueture of the axial divaeaney defeets in 4H- and 6H-SiC in their 
neutral eharge state with S' = 1 eleetron spin. We have previously determined the Dqs parameter 
for hh and kk divaeaneies in 4H Sid^, and now we report it for hh, kiki and ^ 2^2 divaeaneies in 
6 H-SiC. We list the results in Table|^ The results imply that QEl, QE2, and QE6 ODMR signals 
are assoeiated with fci/ci, hh, and ^ 2^2 divaeaneies, respeetively. 
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TABLE I. The calculated hyperfine tensors for the NV center in diamond in the ground (GS) and excited 
(ES) states. The hyperfine constants (A^x, Ayy, A^z) are shown as well as the direction cosine of the largest 
Azz hyperfine consfanf represenfed by angle 6, which is Ihe angle befween fhe direcfion of Azz and fhe 
symmefry axis. The Az hyperfine consfanf is fhe projecfed hyperfine fensor onto fhe symmefry axis. The 
sites are defined by Smelfzer et al^. 


NucleuSsite, state 

Axx (MHz) 

Ayy (MHz) 

A,, (MHz) 

A. 

en 

GS 

3.9 

3.9 

3.4 

3.8 

0 

ES 

-38.5 

-38.5 

-58.1 

-46.0 

0 

l^Ca, GS 

114.0 

114.1 

198.4 

147.6 

71.7 

l^Ca, ES 

44.8 

45.0 

117.5 

77.0 

69.1 

i^Ca, GS 

12.7 

12.8 

18.5 

14.9 

72.0 

^^Ca, ES 

9.6 

9.7 

15.1 

11.8 

73.2 

^^Cb, GS 

11.5 

11.6 

17.0 

13.6 

68.3 

^^Cb,ES 

9.7 

9.7 

15.3 

11.8 

64.8 

^^Cc, GS 

-10.3 

-10.5 

-8.4 

-10.0 

28.1 

^^Cc, ES 

-6.9 

-lA 

-3.6 

-6.2 

13.5 

i^Cd, GS 

-6.8 

-7.2 

-3.8 

-6.1 

71.8 

i^Cd, ES 

-lA 

-7.8 

-5.3 

-6.9 

79.1 

i^Ce, GS 

2.9 

3.0 

4.8 

3.6 

29.3 

^^Ce, ES 

0.7 

1.4 

2.5 

1.6 

23.1 

i^Ce, GS 

4.5 

4.9 

2.9 

4.2 

55.1 

i^Cp, ES 

3.2 

3.8 

4.8 

4.0 

40.1 

^^Cg, GS 

2.1 

2.2 

3.5 

2.6 

75.3 

i^Cg, ES 

1.3 

1.3 

2.4 

1.7 

77.4 

^^Ch, GS 

1.0 

1.0 

2.0 

1.4 

15.0 

i^Ch, ES 

2.0 

2.0 

3.5 

2.6 

13.9 


Proximate nuclear spins of divacancies in GS have been detected by means of ODMRI^, 
labeled as Sina and Sinb by following the labels in Ref.[52l(see Fig.[^. The hh and ^ 2^2 divacancies 
show similar GS hyperfine constants, whereas the kiki divacancy exhibits larger values in 6H SiC 
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TABLE II. A summary of the spin-transition energies for the c-axis-oriented PL6 defect and neutral diva¬ 
cancies in 4H- and 6H-SiC. The parameters are all 20 K parameters, except for the L>es of PL6, where the 
room-temperature value is given. Both Dqs and Des are positive. Comparing the experimental Dqs with 
the calculated Dos (calculated at T=0 K, using the method in Refs. BTl and [60l) allows each spin resonance 
transition in 6H-SiC to be corresponded with its form of neutral divacancy. 


Defect 

Dos (GHz) 

(GHz) 

Des (GHz) 

4H: hh 

1.336 

1.358 

0.84 

4H: kk 

1.305 

1.320 

0.78 

4H: PL6 

1.365 

- 

0.94 

6H: hh 

1.334 

1.350 

0.85 

6H: kiki 

1.300 

1.300 

0.75 

6H: ^2^2 

1.347 

1.380 

0.95 


(see Table [nil). The kiki in 6H-SiC and kk in 4H-SiC as well as hh configurations in 6H- and 4H- 
SiC show similar values. These trends support the assignment of QL1,2,6 based on the calculated 
Dos parameters. 

Finally, we calculated the corresponding hyperfine tensors in the ES (see Table jl^. Within the 
accuracy of measurements, the ES shows signatures of C^v symmetry. However, a full character¬ 
ization of the divacancy’s ES is beyond the scope of this study. We approximate the ES by the 
constraint DET procedure and fix the C^v symmetry in our calculations to obtain the ES hyperfine 
tensors. 


V, RESULTS AND DISCUSSION ON DYNAMIC SPIN POLARIZATION PROCESSES 

In this section, we present our results and conclusions on the DNP process for the NV center 
in diamond and for the divacancy in 6H-SiC. Eirst, we start with the study of the NV center in 
diamond, which is the most thoroughly investigated defect system exhibiting DNP. We show that 
the new model is capable of reproducing both theoretical and experimental curves. Nevertheless, 
it provides a new and deeper insight to the features of the DNP mechanism and to the physics of 
the defect. After this, we apply our model for different configurations of the divacancy in 6H- 
SiC. Our results reveal the importance of the electron spin coherence time in the excited state in 


21 







TABLE III. The calculated hyperfine tensors for nuclear spins that are proximate to divacancies in 6H-SiC 
in the ground (GS) and exited (ES) state. The hyperfine constants {A^x, Ayy, A^z) are shown, and also 
the direction cosine 0, which is the angle between the direction of Azz and the symmetry axis. The Az 
hyperfine consfanf is fhe projecfed hyperfine fensor onfo fhe symmefry axis. The afom labels are shown in 
Eig. Since PE6 in 4H-SiC has nol yef been identified, we applied fhe calculated hyperfine fensors in fhe 
^ 2^2 divacancy configurafion in fhe DNP simulations. In our experience, fhe difference in fhe measured 
or calculafed hyperfine consfanfs are wifhin I MHz, and a I MHz inaccuracy does nof alfer fhe simulation 
resulfs. 


Nucleus 

Site 

Conf. 

St. 

A,, [MHz] 

Ayy [MHz] 

A,, [MHz] 


A^ [MHz] 

29Si 

Siiib 

6H: hh 

GS 

9.8 

8.6 

10.7 

69.5 

9.6 

29Si 

Siiib 

6H: hh 

ES 

9.8 

9.3 

10.4 

63.0 

9.8 

29Si 

Siiib 

6H: kiki 

GS 

10.7 

9.8 

11.5 

70.8 

10.5 

29Si 

Siiib 

6H: kiki 

ES 

10.2 

9.5 

10.8 

60.7 

10.1 

29Si 

Siiib 

6H: k2k2 

GS 

9.9 

8.7 

10.8 

69.4 

9.7 

29Si 

Siiib 

6H: fc2/c2 

ES 

10.4 

9.7 

11.0 

64.1 

10.3 


understanding DNP. 


A. NV-center in diamond and a single adjacent nuclear spin 

Previous modelP^ sueeessfully explained experimental results on the dynamie nuelear spin 
polarization of a nuelear spin of the NV-eenter in diamond. In this speeial case, the hyperfine 
tensor of the excited and ground states are symmetric (Table |^. Therefore, when the magnetic field 
is well aligned with the axis of the defect the spin Hamiltonian becomes rather simpleP^^. Here, 
we reassess these experimental and theoretical results in the framework of our extended method. 

To calculate the magnetic field dependence of the nuclear spin polarization, we chose param¬ 
eters for our model as follows. The parameters of the hyperfine tensor of the ground and excited 
state are determined by our first principles calculations. Table |Ij The D parameters for the ground 
and excited state zero-field-splitting tensors were set to 2.87 and 1.42 GHz, respectively, in accor¬ 
dance with experiment. For the excited state the lifetime and the rate of non-radiative decay, we 
used the experimental values t^s = 12 ns and T = 0.3. The jj, parameter in our model determines 
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FIG. 4. A comparison of calculated and measured dynamic nuclear spin polarization P of a nucleus 
of the NV-center in diamond as a function of the external magnetic field B. The result of our calculation 
is depicted with a thick red solid line, while the previous theoretical and experimental resultJ^ are de¬ 
picted with thin black solid line and black points, respectively. The calculated ESLAC resonance peak, at 
^ESLAC = 516 G, resembles the reported oned^. However, in our calculation the ground state processes 
are also taken into account, producing an additional peak at Pgslac = 1020 G. 

the features of the GSLAC resonanee peak. Given the laek of experimental data for this peak, we 
set /i = 0.25. The shape of the ESLAC polarization curve is determined by two parameters in 
our model. Parameters k and u take into account the effect of the nuclear spin-lattice relaxation 
and the excited state’s electron spin decoherence, respectively. These parameters were fit to the 
experimental curves P{B) and P{Qb)- The optimal values are u = 0.13 and k = 3.15 x 10“^. 

In order to indicate the sensitivity of the fitted curve to the variation of the free parameters, we 
depicted curves of modified v and k parameters in Fig. Q. As can be seen the theoretical curve 
moderately but non-linearly varies with the changes of the fitting parameters. 

We note that the fitting parameter k has a smaller value in our model than parameter = 
0.0027 of the previous model with same definition. The difference is due to the effect of the 
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FIG. 5. Parameter dependence of the theoretical nuclear spin polarization curve P{B) for the case NV 
center in diamond including an nuclei, (a) and (b) show u and k parameter dependence of the theoretical 
curve, respectively. In both cases, the red thick curve corresponds to the optimal parameter setting, see 
Fig. @, while the edges of the light purple filled area correspond to ±60% change of the parameters. 
Variation of parameter ^ has similar effect on the GSLAC peak as u has on the ESLAC peak. 

ES electron spin decoherence. In the previous model, this effect is not included, while in our 
model it is taken into account explicitly. The electron spin decoherence shortens the effective ES 
lifetime, thereby reducing the probabilities of nuclear spin flips and the polarization. This finding 
immediately draws attention to the importance of the electron spin decoherence. 

Our extended model reproduces well the reported magnetic field dependence of the nuclear 
spin polarizatioiJ^ at the vicinity of I^eslacj see Eig.|^ The position of the maximum polarization 
is at i? = 516 Gauss. It is slightly shifted due to the different hyperfine tensors utilized in the 
two calculations. On the other hand, in our model, the ground state processes are simultaneously 
described that produces an additional resonance peak at I^gslac = 1020 Gauss. As the ground 
state hyperfine interaction is weaker than the excited state interaction, the GSEAC resonance peak 
is much sharper that at the ESEAC. 

We also calculated the magnetic field angle dependence of the degree of P, see Eig. which 
was not considered by previous models. The theoretical curve agrees well with the result of the 
experimental measurement^. The polarization decays very quickly as the angle of the magnetic 
field increases. Our model makes it possible to investigate these observations thoroughly. Eirst 
of all, as the symmetry of the hyperfine tensor is reduced by the non-zero angle of the magnetic 
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(a) 



FIG. 6 . The magnetic field angle dependence of the nuclear spin polarization, (a) shows the calculated 
(solid line) and measured^ (points) spin polarization of a nucleus of an NV-center in diamond as a 
function of the angle of the external magnetic field and the 6 * 3 ^ axis of the defect at i? = 516 Gauss. 
Our calculation accurately reproduces the experimental observations, (b) The angle of the nuclear spin 
is depicted as function of the angle of the magnetic field. As the misaligned magnetic field reduces the 
symmetry, the stationary state of the nuclear spin points out of the axis. 

field and the symmetry axis of the defeet in the stationary state, the nuelear spin is not parallel 
with either the symmetry axis of the defeet or the magnetie field. The angle Op of the preferential 
direetion of the polarization linearly depends on the angle Op of the magnetie field, see Fig.|^ The 
ratio of the two angles is Op/Op = 4.03. 

In Fig. 1^ we depieted the magnetie field dependenee of the nuelear spin polarization for the 
ease of a misaligned magnetie field. The angle of the magnetie field and the symmetry axis of the 
defeet is set to 1°. The ESLAC resonance peak is strongly reduced, whereas the GSLAC resonance 
peak almost disappears. 

Next, we investigate the role of different spin rotation processes when the symmetry is reduced. 
In Fig.[^ we compare the time dependence of the ES spin flipping processes for the case of aligned 
and misaligned magnetic fields. In the former case, only one spin rotation mechanism can be 
observed, |0 J,) —?• |—1 in accordance with the previous modelP^^. By tilting the direction of 
the magnetic field, the maximal probability of this process decreases while other flipping processes 
appear simultaneously with high probability (see Eig. [^b) ). All of these mechanisms reduce the 
maximal polarizability of the nuclear spin. Eor instance, |0 t) —t | — 1 i) represents a driving force 
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FIG. 7. The calculated degree of dynamic nuclear spin polarization of the nucleus of the NV-center in 
diamond for the case of misaligned magnetic field. The angle of the magnetic field and fhe axis of fhe 
defecl is sef fo 1°. (a) The polarizafion of fhe nuclear spin as a function of fhe sfrengfh of fhe magnefic field. 
If shows a sfrong reducfion in fhe polarizabilify compared wifh Fig. due fo fhe misalignmenf of magnefic 
field. Interestingly, fhe ground-sfafe dynamic nuclear spin polarization almosf completely disappears, (b) 
The angle of fhe nuclear spin and fhe axis of fhe defecl as a funclion of fhe magnetic field. 


towards nuclear spin polarization in the opposite direetion, while |0 ],) —)■ |0 t) represents the 
preeession of the nuclear spin that reduees the polarization regardless the direction of the nuelear 
spin. The interplay of these effeets is responsible for the reduction of the ES polarizability. 

From Fig. it is elear that the ES lifetime and the length of intact evolution time play an im¬ 
portant role in DNR By fitting the theoretieal eurve to the experimental data P(0 b), we obtain 
the average length of the intaet evolution time, whieh turns out to be only 14% of the ES lifetime. 
This short time suggests substantial eleetron spin deeoherenee in the excited state. In sueh eircum- 
stanees, the spin rotation proeesses have a longer timeseale than the net evolution time. There are 
two main eonsequenees of this. First, the nuelear spin flipping probabilities are largely redueed, 
sinee the spin rotations are suppressed. Seeond, there is a preferenee over fast proeesses. For 
example, in the case of misaligned magnetic field, the fastest meehanism, exhibiting rapid oscil¬ 
lations, can flip the spins with the highest probability during the short time evolution. Thus, at 
the ESFAC |0 ],) —)■ |—1 t) proeess is still dominant eausing non-zero nuelear spin polarization. 
On the other hand, by elongating the ES evolution time, this preferenee relaxes and other slower 
proeesses ean have more pronouneed effeets (see Fig.[8]^b) ), whieh would greatly ehange the nu- 
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FIG. 8. The time dependence of spin flipping probabilities in the excited state of the NV center in diamond 
with a nucleus, (a) The magnetic field is well aligned with the axis of the defect. In this case, only one 
spin rotation can occur, |0 J,) —)• | —1 t)^ represented with red (thick gray) line. The black and blue lines 
with light gray and blue filled areas represent the exponential decay of the excited state lifetime and the 
effective evolution time, respectively (see text for more explanation). The mean evolution time can be seen 
to be much shorter than the periodicity of the oscillatory probability, (b) The magnetic field is misaligned 
by 1°. In this case, different kinds of spin rotations occur that lower the polarizability of the nuclear spin 
(Not all processes are depicted). 

clear spin polarization pattern. This ean be the reason for the disappearanee of the nuelear spin 
polarization at the GSLAC (see Fig. |^), where the ground state’s evolution time is mueh longer 
than that of the exeited state. 

Finally, we emphasize the strong relation between the magnetie field angle dependenee of the 
nuelear spin polarization and the eleetron spin deeoherenee in the ES. Shorter eoherenee time 
reduees the maximal polarizability but the deeay of the polarization with respeet to angle ©b of 
the magnetie field is slower. 


B, NV-center in diamond and first neighbor nuclear spin 

In this seetion, we study the dynamie nuelear spin polarization of a nuelei in the first 
neighbor shell, Cq site, of the NV eenter, where both the ESLAC and GSEAC polarization was 
observecP^*^. In the latter ease, partieular features appear due to the non-symmetrie hyperfine 
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FIG. 9. Calculated magnetic field dependence of the dynamic nuclear spin polarization of a single nucleus 
in the first neighbor shell of the NV-center in diamond, (a) The polarization of the nuclear spin as a 
function of the strength of the magnetic field, (b) The angle of the nuclear spin and the C^y axis of the 
defect as a function of the magnetic field. 


interaction, which show a potential for new applications in the area of sensitivity-enhanced nuclear 
magnetic resonance and spintronic^. 

As a N atom with non-zero nuclear spin is an inherent part of the NV center, the inclusion of 
an adjacent nucleus necessarily results in a three-spin system. Due to the small gyromagnetic 
ratio of the nuclei, the direct nuclear spin-nuclear spin interaction is negligible. Substantial in¬ 
teraction is, however, mediated by the electron. In this case, the hyperfine interaction with the 
nitrogen atom appears as an additional magnetic field dependent electron spin relaxation effect 
for the two-spin system of the nucleus and the electron. On the other hand, in the stationary 
state of the system, the nitrogen nuclear spin is highly polarized, while the populated spin state of 
the nitrogen nucleus is only weakly coupled to the electron spin. Therefore, we neglect the effect 
of the nitrogen nuclear spin on the polarizability of the nucleus. We investigate the system 
as a two-spin system within our model. Description of three-spin DNP systems can be found in 
Ref. HU 

In the calculation, we used the hyperfine parameters presented in TableOther parameters are 
the same as in the previous section. The reason for keeping the parameters is that the measurements 
on the ESLAC polarizability of and nuclei were carried out on same sample!^. 

As can be seen in Fig.|^ the predicted maximal nuclear spin polarization of 70% at ESLAC is 
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slightly below the experimentally observed Since the measurements were carried out on 

single centers, the local environment is not necessarily the same, which may be the reason for the 
differences. By assuming weaker spin dephasing effects, higher polarization can be achieved in 
the calculations. 

The oscillatory features of the calculated nuclear spin polarization of the nucleus at GSLAC 
(see Fig. |^a) ) reproduces well the experimental observations by Wang et alW^. In our case the 
nuclear spin polarization reaches 83% at 950 Gauss before it rapidly drops at GSLAC. Right after 
the GSLAC, the polarization reverses and reaches -11% at 1050 Gauss. It then increases again 
and reaches 62% at 1170 Gauss. This feature has been thoroughly investigated and understood 
by Wang et alW^. In short, due to non-symmetric hyperfine interactions, unusual combinations 
of ladder operators appear in the spin Hamiltonian (see Appendix). From Fig. [^a) and (b) it is 
clear that the reverse of the polarization is due to the dominance of process |0 t) | —1 i) over 
the most common process |0 ),) —)■ | —11) in the ground state. This dominant process occurring 
only at well-defined external magnetic fields is governed by the 5_J_ operator combination that 
resonantly mixes state |0 t) and | — 1 ),) at a certain magnetic field slightly higher than I^gslac- The 
resonance appears as a dip due to the inverse polarization. Since the hyperfine constant correspond 
to term J_ is smaller than the one correspond to the regular term S-i+, the dip is sharper than 
the resonance peak due to the term J+. In summary, GSLAC polarization curve built up from a 
wide resonance peak, due to term ^_/+, and a slightly shifted sharper dip, due to term J_. 

By looking at the polarization near the ESLAC (Fig.|^, one may recognize similar oscillation 
of the polarization. However, here the pattern is not as developed as for the GSLAC resonance. 
From the discussion presented in the previous section, it is clear that the short ES lifetime and the 
fast electron spin dephasing effect cause the difference. The hyperfine constant of term is 
larger than that of which implies slower spin rotation processes for the latter term. As the 

short evolution time reduces the flipping capability of slow processes, the dip corresponding to the 
slower S-I- process is effectively suppressed. By manually increasing the ES lifetime, we can 
see a more developed dip. 

Einally, the magnetic field dependence of the polarization direction of the nucleus exhibits 
similar pattern at the GSEAC (see Eig. |^(b) ) as for the case of misaligned magnetic field. 
One can argue that the two types of symmetry reduction of the spin Hamiltonian result in quali¬ 
tatively the same phenomena. This realization may provide a way to engineer the magnetic field 
dependence of the polarization of non-symmetrically placed nuclei. 
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C. NV-center in diamond and remote nuclear spins 


The study of the polarizability of remote nuelear spins provides an opportunity to have a 
deep insight into the parameter dependence of the maximal achievable nuclear spin polarization in 
DNP processes, which is thoroughly examined in the first subsection. Our conclusions are tested 
on the available experimental results in the second subsection. 


1. Parameter dependence of nuclear spin polarization 

From the applications point of view, it is important to understand the degree to which the spin 
polarization of the electrons can be transferred to the nuclei by DNP and what the limiting fac¬ 
tors are. Remote nuclear spins couple weakly to the electron spin of the defect, and thus their 
polarizability is assumed to strongly depend on the details of the hyperfine tensor. In this section, 
we aim to understand the role of different hyperfine parameters in the determination of the maxi¬ 
mal achievable nuclear spin polarization of more remote nuclei that couples relatively weakly to 
the electron spin of an NV center. Previous publications have already addressed this issue to an 
extenP®^. However, a satisfactory conclusion has not been achieved so far. The most commonly 
mentioned reasons for the varying polarizability of remote nuclei are the varying angle of the 
principal axis of the ES hyperfine tensor and the symmetry axis of the defect as well as the varying 
angle of the principal axis of the ES and GS hyperfine tensor. On the other hand, as we will show, 
the strength of the coupling may have the most important role. 

In this section, we investigate a set of (/ = 1/2) nuclear spins interacting with the electron 
spin of the NV-center in diamond at ESEAC condition. We examined the achievable nuclear 
spin polarization as a function of the strength of the coupling, the anisotropy of the hyperfine 
interaction, the direction cosines of the ES hyperfine tensor, and the difference of the direction 
cosines of the ES and GS hyperfine tensors. In the calculations, we arbitrarily vary = A^y = 
and Aff = in the range of 1-24 MHz and 6^^ and 6^^ in the range of 0-90° to map the 
parameter space, while we kept fixed other parameters of the model. The free parameters were 
initialized as k = 0.0015, u = 0.5, and p = 1. 

Eirst, we investigate the effect of GS processes on the ESEAC polarizability of weakly coupled 
nuclear spins. The GS-spin Hamiltonian plays an important role, since the system spends most 
of the time in GS. In our model, this Hamiltonian determines the direction of nuclear spin polar- 
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ization. This suggests that the details of the GS hyperfine interaetion can affect the polarizability. 
Our calculations, however, show that the nuclear spin polarization direction is parallel with the 
axis of the defect, i.e. the nuclear spin polarized parallel to the electron spin, regardless the details 
of the GS hyperfine tensor. Thus, we can say that the difference of the principal axis of the GS 
and ES hyperfine tensor, which is thought to play an important role, as well as other parameters 
of the GS hyperfine tensor, do not affect the ESLAC polarizability of the spin of remote nuclei. 
One may prove this observation based on simple analytical arguments. Consider the usual basis 
|0 t), |0 i), |—1 t)> and | —1 i), where the quantization axis of the electron and nuclear spins are 
parallel to the symmetry axis of the defect. Ear from the GSEAC, the large zero field splitting 
of the electron spin states suppresses the mixing of the states, while at the GSEAC the mixing is 
more enhanced due to the GS hyperfine coupling. On the other hand, for the case of weak GS 
interactions, the mixing and other effects of the interaction are simply negligible at the ESEAC, 
since A ^ 0(10 MHz) C 1400 MHz Dos — J^es- Thus, the state |0 t) and |0 |) are eigenstates 
of the GS spin Hamiltonian within good approximation. Their long GS lifetime makes them good 
target states for DNP polarization. 




FIG. 10. Parameter dependence of the dynamic nuclear spin polarization of nuclei adjacent to a NV 
center in diamond, (a) The polarization is plotted against the angle of Csv axis and the principal axis of the 
ES hyperfine tensor for different values of coupling strength A±. (b) The polarization is plotted against the 
coupling strength A±. The hyperfine parameters are given in MHz unit and refer to ES interactions. 

To investigate the dependence of the polarizability on the details of the ES hyperfine tensor, 
we depict the maximal polarization of nuclear spin as a function of the principal axis of 
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the hyperfine tensor for different eoupling strength A± in Fig. [T^ (a). As one ean observe, the 
angular dependenee of the polarization is not substantial in most oases. However, it depends on 
the anisotropy of the hyperfine interaetion. For example, when the interaetion is isotropie, i.e. 

= A\\, no angular dependenee is olearly observed. However, when the interaetion is highly 
anisotropio more pronounoed angular dependenee ean be observed, see the ease of A^/A\\ = 1/8 
or A^/A\\ = 3. To understand the observed features, we need to take into aeoount three effeots: 
First, the value of the off-diagonal elements, whioh are responsible for the spin flipping prooesses, 
are funotions of angle 9 as 



(24) 


For instanee, the mixing term inereases with angle 9 as A\\ > A^, whioh results in inereasing 
polarizability for small 9, see the ease of A± = 1 MHz and 2 MHz in Fig. [T^(a). The opposite 
effeot ean be observed for Ay < A±. This effeot inereases with the anisotropy of the hyperfine 
tensor. Seeond, as the angle 9 inereases, the symmetry of the spin Hamiltonian beeomes more 
and more reduoed, thereby strengthening the nuolear spin depolarization proeesses (see Fig.[^(b) 
). Third, the effeot of these prooesses depends on the net evolution limited by the lifetime of the 
exoited state and the eleotron spin ooherenoe time in the exoited state. The interplay of the above 
desoribed two basio meohanisms determines the hyperfine angle dependenee of the polarizability. 

Next, we investigate DNP as a funetion of the strength of the hyperfine eoupling. We depiet 
the maximal polarization against parameter A± in Fig. [I^b) with fixed A\\/A±_ ratio and angle 
9. As one ean see, the polarization quiekly drops as the hyperfine interaetion deereases below 
10 MHz. The reason of the reduotion of the polarizability is the finite evolution time in the ES of 
the dynamieal eyele. As the strength of the eoupling deereases, the rate of the spin flipping proeess 
deereases as well. Due to the finite lifetime, the time is not suffieient to flip the spins. On the 
other hand, the eharaeter of the eurve depieted on Fig. [T^b) strongly depends on the length of the 
effeetive evolution time and thus on the strength of the eleetron spin dephasing effeets. We ean say 
that the finite lifetime and the short eleotron spin ooherenoe time of the ES limits DNP of remote 
nuolear spins. Shorter evolution times oauses the drop of the polarizability at larger hyperfine 
eouplings. Here, we mention that the enriehment of the nuolear spin oonoentration effeotively 
deereases the eleetron spin eoherenee time, thus suppressing DNP prooesses of remote and nearby 
nuolear spins. This meehanism ean have important eonsequenees on the NMR applieations of 
these DNP prooesses, and may suggest the importanoe of implementing dynamieal deooupling 
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protocols designed to enhanee DNR 

To summarize our findings on the hyperfine parameter dependenee of the DNP proeess we ean 
say that, at a given external eondition, the polarizability of remote nuelear spins is determined 
by their A± parameter in the first order. A diverse seeondary eontribution is governed by the 
relationship of the three parameters A±, Ay, and 0. The details of the GS hyperfine interaetion 
do not affeet the polarizability. And finally, the eleetron spin eoherenee time limits how far the 
eleetron spin polarization ean be transferred by ESLAC DNP meehanism. 

2. Comparison of theory and experiment for the case of remote nuclei 

In this subseetion, we eompare the results of our ealeulations with the available experimental 
measurements for the ease of distant nuelei around an NV eenter in diamond (see Fig. [^. In 
our ealeulations, we used the ealeulated hyperfine eoupling parameters listed in Table while the 
zero-field-splitting parameters, the ES lifetime, and the non-radiative deeay rate are the same as in 
the previous seetions. The free parameters were speeified as k = 0.003, p = 0.9, and z/ = 1. In 
the ealeulations the magnetie field was set to B = 510 Gauss as given in the experiments. 

The ealeulated polarization is in good agreement with the experimental measurements by 
Smeltzer et al^. However, the agreement is worse for the ones obtained by Dreau et al^ (see 
Fig.[^. One should keep in mind that there are many faetors that may limit the direet eomparison 
between theory and experiments. Beside the experimental uneertainties, the ealeulated ES hyper¬ 
fine parameters have systematie uneertainties due to the approximate theory and the diffieulties of 
the deseription the ES-eleetron eonfiguration. 

By looking at the theoretieal result in Fig. as well as at the hyperfine parameters in Table 
the basie roles for the polarizability, defined in the previous subseetion, ean be reeognized. How¬ 
ever, there is an additional faetor one should take into aeeount. Different sites exhibit different 
hyperfine eoupling to the eleetron spin of the NV eenter. By applying the same magnetic field for 
all the different sites the resonance condition is not satisfied equally for all the nuclei spins. As the 
hyperfine eoupling is relatively weak and the resonanee peaks are sharp, the polarization pattern 
among the different sites ean sharply depend on the value of the external magnetie field. This may 
be a reason for the substantial deviation between the two available experimental measurements. In 
our ealeulations, we found that the variation of A5 = ±5 Gauss ean result in 5-25% variation of 
the observed polarization among the different sites. Our previously defined rules apply to maxi- 
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FIG. 11. The degree of dynamic nuclear spin polarization of nuclei at various sites around a NV 
center in diamond. The calculated polarizations (the red line with stars) are compared with the available 
experimental results. Blue squares and green triangles show the experimental results by Smeltzer et a/.l^and 
Dreau et alW^, respectively. The sign of the experimental polarization is changed for sites C and D, see text 
for more explanation. Sites E and F and sites G and H could not be separated in the experiment due to their 
similar GS hyperfine parameters. Thus, the measured polarization can be the mixture of the polarization of 
the two sites. As sites F and H with three-fold degeneracy are more scarce than the corresponding E and 
G sites with six-fold degeneracy, F and H data points are represented with empty symbols. Because only a 
few centers with nuclei on a specific site were investigated, it is possible that only the most probable 
sites, E and G, were measured. 

mally achievable polarization and thus to perfeet resonanee eonditions. The imperfeet resonanee 
eondition overshadows the trends dietated by our findings. 

Our ealeulations show positive polarization for both positive and negative hyperfine parameters. 
However, in previous manuseripts, negative hyperfine parameters were assoeiated with negative 
polarizatiorP^*^. We think that the diserepaney eomes from a misinterpretation of the sign of the 
polarization in the experiments. For positive hyperfine parameters, the population of the lower 
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energy eigenstate represents positive nuelear spin polarization. On the other hand, for negative 
hyperfine parameters, positive polarization eorresponds to the population of the higher energy 
eigenstates. In the measurements, the population of the higher energy eigenstate is interpreted 
as negative polarization. However, it is actually due to the change of the sign of the hyperfine 
parameter. Thus, in Fig.[^ we depiet the experimental polarization of C and D sites with positive 
sign, in agreement with our simulations. 


D. Divacancy in 6H-SiC and adjacent ^®Si nuclear spin 

Recently, dynamic nuclear spin polarization of the hh and kiki configurations of the divacancy 
in 6H-SiC and the related PL6 center in 4H-SiC were observed^. In this section, we use our model 
to reproduce the measured magnetic field dependence of the nuclear spin polarization at the vicin¬ 
ity of the ESLAC and to use the fitting parameters to gain important insight to the DNP processes. 
Here, we restrict ourselves to study of ESLAC DNP of the different divacancy configurations. 

Considering the experimental result^ and the calculated hyperfine parameters of the different 
configurations, see Table |II^ an interesting observation can be made. Although, the hyperfine 
parameters of these configurations differ only a few MHz, the measured nuclear spin polarization 
vary between 60% and 99% at ESLAC, see Eig. As we will show in this section, such large 
variation of the polarizability can only be explained by the different ES electron spin coherence 
times of the different configurations. This observation further emphasizes the role of these effects 
in the DNP process. 

TABLE IV. The utilized excited state zero-field-splitting parameter Des as well as the measured^ rate of 
non-radiative decay T and excited state lifetime tes are listed for hh and kiki configuration of divacancy 


in 6H-SiC and the related PL6 defect in 4H-SiC. 


Configuration, SiC polytype 

Des [MHz] 

r 

tes [ns] 

hh, 6H 

930 

0.15 

15 

kiki, 6H 

915 

0.16 

15 

PL6, 4H 

955 

0.33 

14 


The hyperfine parameters, used in the DNP calculations, were determined by our ab initio 


calculations, see Section IV Recent experiments revealed sensitive temperature dependence of 
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FIG. 12. The adjusted theoretical curves (solid lines) and the experimental measurement(points) on 
the dynamic nuclear spin polarization of different divacancy configurations and related defects. The panels 
show the QLl, QL2 and PL6 photoluminescence centers (ordered from left to right). The first two centers 
are assigned to the kiki and hh configurations of divacancy in 6H-SiC (see text for further information). 
The unidentified PL6 center is modeled by hyperfine parameters of k 2 k 2 configuration in the theoretical 
simulations. 

the excited state zero-field-splitting of the divacanc)^^. Since the DNP measurements were carried 
out at various temperature^, we also adjusted Des in our DNP simulations. The E^s, parameter is 
set to zero in all cases. The D parameter together with the measured rate of non-radiative decay T 
and ES lifetime Tes are listed in Table The fitting parameters correspond to the different type 
of spin relaxation processes are collected in Table fV] 

The adjusted theoretical curves and the experimental measurements are depicted together in 
Fig. for hh and kiki configuration of divacancy in 6H-SiC and PL6 center in 4H-SiC. Good 
agreement is achieved for all the cases. We would like to emphasize that the theoretical curves 
are adjusted by the variation of parameters that are related to spin decoherence and spin-lattice 
relaxation processes. Through the effect of the variation of these parameters, the effect of different 
spin relaxation processes may be understood. The k parameter, corresponding to the rate of spon- 


36 







TABLE V. Fitting parameters k and v for hh and kiki configuration of divacancies in 6H-SiC and the 
related PL6 defect in 4H-SiC. The latter parameter determines the effective evolution time which is then 
compared with the experimentally measured^ coherence time r|. In addition, the ambient measurement 
temperature T is given. 


Configuration, SiC polytype 

K 

V 

Tes = [ns] 

r| [ns] 

T [°K] 

hh, 6H 

8.25 X 10-5 

0.092 

1.38 

4.0 

100 

kiki, 6H 

8.25 X 10-5 

0.045 

0.68 

1.3 

100 

PL6, 4H 

3.14 X 10-^ 

0.362 

5.1 

5.0 

300 


taneous nuclear spin flipping and thus to the nuelear spin-relaxation time, primarily determines 
the deeay of the polarization eurve away from the ESLAC resonanee. Seeond, it lowers the max¬ 
imal aehievable polarization at resonanee. The v parameter eorresponds to the effeet of eleetron 
deeoherenee in the ES. It mainly determines the height of the resonanee peak, thus the maximal 
polarization at ESEAC. 

Spin relaxation proeesses are taken into aeeount rather approximately in our model, whieh 
might deseribe either the intrinsie properties of the defeets or the different environment of the 
defeets. Considering the fitting parameters in Table [V] one ean see that k has larger value for PE6 
in 4H-SiC than for hh and kiki in 6H-SiC, suggesting faster nuelear spin-lattiee relaxation for 
PE6. The reason ean be threefold: (i) the measurements on the PE6 eenter were earned out at 
room temperature while in the other eases at 100 K; (ii) the 4H sample is differently proeessed 
than 6H sample with possibly ereating different number of disturbing spins the erystals; (iii) PE6 
is preferentially loeated relatively elose to the surfaeeP^ where the eoneentration of non-eoherent 
eleetron spins may be larger. In Table |V| one also finds the values of the fitting parameter u, 
together with the estimated effeetive evolution time with the trend of kiki < hh < PE6. This 
trend is good agreement with the one that ean be observed from the measured eleetron dephasing 
time T 2 , see Table fV] governing the length of the effeetive evolution time in the exeited state 
in our model. This result indieates that the aehievable high (low) ESEAC polarizability for PE6 
(kiki) is due to the longer (shorter) eleetron eoherenee time in the ES. 
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VI, SUMMARY 


We developed a model that provides a detailed description of the dynamic nuclear spin polariza¬ 
tion process of point defects with a high ground-state electronic spin and adjacent nuclear spins. It 
takes into account many microscopic features and processes that previous models have overlooked 
or not integrated: the effect of the electron decoherence in the ground and excited state, the short 
lifetime of the excited state, the direction of the nuclear spin polarization, the anisotropy of hy- 
perfine tensors, misaligned magnetic fields, and the simultaneous flipping of ground and excited 
state’s spin flipping processes. We applied our model for different cases of the NV center in dia¬ 
mond and divacancies and related defects in 4H- and 6H-SiC, validating our model and showing 
its generality. Additionally, in many cases, the model provided a new insight to the microscopic 
mechanisms behind the examined phenomena and allowed us to draw important conclusions. We 
show that the electron decoherence in the ES plays an important role in DNP process, through the 
shortening the ES evolution time. It can limit the maximal achievable ES nuclear spin polariza¬ 
tion, as we showed for the divacancy in SiC. By promoting fast nuclear spin rotation processes 
and suppressing slower ones, the electron decoherence protects the ESEAC resonance from sec¬ 
ondary spin rotation mechanisms that appear for anisotropic hyperfine interactions. On the other 
hand, similar secondary processes produce more complicated fine structure at GSEAC, due to a 
longer evolution time permitted by the longer lifetime of the GS in the optical cycle. By studying 
the polarizability of remote nuclear spins, we addressed the long-standing question concerning 
the hyperfine parameter dependence of the maximal achievable polarization. As we showed, the 
strength of A± primarily governs the polarizability of remote nuclei, while the ES electron spin 
decoherence is an important factor as well. Throughout our investigations, electron spin decoher¬ 
ence appeared as an important, sometimes major, limitation of the polarizability of the nuclear 
spins. Therefore, any treatment that affects the point defect electron spin coherence times, like 
isotope enrichment, can drastically affect the DNP process as well. 
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Appendix: Spin Hamiltonian matrices 


Here, we specify the spin Hamiltonian matrices for the coupled system of a S' = 1 electron 
spin and an / = 1/2 nuclear spin. The spin Hamiltonian consists of four terms (see Eq. (1) in the 
main text). The corresponding four matrices, as written in the basis of |-|-1 t), |-l-l i), |0 t), |0 ),), 
| — 11 ), and | —1 i) states, where the spin quantization direction 2 : is parallel to the G 3 axis of the 
defect, are listed below. The occurrence of non-zero off-diagonal elements, which are responsible 
for spin rotation processes among the basis states, are discussed in each cases. 


• The zero-field interaction term 


H.fs{D,E) 


/d 0 0 0 E o\ 
0 D 0 0 0 E 
0 0 0 0 0 0 
0 0 0 0 0 0 
E 0 0 0 D 0 
\0 E 0 0 0 Dj 


(A.l) 


where D and E are the zero-field-splitting parameters. For point defects of symmetry. 


the E parameter is zero. Thus there are no off-diagonal elements in the spin Hamiltonian 
due to the zero-field interaction. 


Electron spin Zeeman term 
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(A.2) 
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where B is the strength of the external magnetic field, is angle of the C 3 axis of the defect 
and the external magnetic field, ge is the g-factor of the electron spin, and /xb is the Bohr 
magneton. 

As can been seen, off-diagonal matrix elements appear for 6^ 7 ^ 0. These elements 
correspond to the electron spin ladder operators S+ and S- and couple \Ms,Mi) and 
\Ms ± I,Mi) states. 


Nuclear spin Zeeman term 


H'sziB, 9 b) — 


(cosOb sin 05 0 
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cos 9b 0 
0 cos 9 b 

0 sin 6*5 
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0 
0 
0 
0 

sin 9 b 
- cos 6*5y 


(A.3) 


where B is the strength of the external magnetic field, 6*5 is angle of the C 3 axis of the defect 
and the external magnetic field, is the g-factor of the nuclear spin, and /tn is the nuclear 
magneton. 


Similarly to the electron spin Zeeman term, off-diagonal matrix elements are present for 
6 *b 7 ^ 0. Here, these off-diagonal elements correspond to the nuclear spin ladder operators 
/+ and /_ and couple \Ms,Mj) and \Ms,Mj ± 1) states. 


• Hyperfine interaction term 


Bhyp(^Axxj -‘^yyi Azzj 9) 
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a = Azz cos^ 9 + Axx sin^ 9 
b = (Azz — ^xx) cos 9 sin 9 
c± = Axx cos^ 9 -f Azz sin^ 9 ± Ayy, 


(A.4) 


(A.5) 

(A. 6 ) 

(A.7) 
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where A^x, Ayy, and define the hyperfine coupling strength and 6 defines the polar angle 
of the principal axis of the hyperfine tensor. 

For the case of 6 * = 0 and Axx = Ayy = A±, the non-zero off diagonal matrix ele- 
ments ^c_|_ = correspond to the ^±/=p operator combinations. These elements 

mix \Ms,Mj) and \Ms ± l,Mi =f 1) states. When Axx ^ Ayy but 6 * = 0, new off-diagonal 
elements are present in connection to the S±I± operators. The coupling strength is ^c_ = 
^ {Axx — Ayy) and the coupled states are \Ms,Mi) and \Ms ± l,Mi ± 1). When 0 7 ^ 0 
and 9 ^ where n is an integer, precession like spin rotation process can appear due to the 
S±Iz and Sj± operators that couple state similar to the electron and nuclear spin Zeeman 
effects, respectively. 
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